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> : ABSTRACT 

^ I A thorough study is performed of the analytical properties of the fermion determi- 

t^^ I nant for the case that the time components of (axial) vector fields do not vanish. For 

this purpose the non-Hermitian Euclidean Dirac Hamiltonian is generalized to the whole 



0\ I complex plane. The Laurent series are proven to reduce to Taylor series for the corre- 

1^ ' sponding eigenvalues and -functions as long as field configurations are assumed for which 

I I level crossings do not occur. The condition that no level crossings appears determines the 

qj I radius convergence. However, the need for regularization prohibits the derivation of an 

'^ ' analytic energy functional because real and imaginary parts of the eigenvalues are treated 

differently. Consistency conditions for a Minkowski energy functional are extracted from 
r> ' global gauge invariance and the current field identity for the baryon current. Various 

d '_ treatments of the Nambu-Jona-Lasinio soliton are examined with respect to these con- 

ditions. 

Motivated by the studies of the Laurent series for the energy functional the Euclidean 
action is expanded in terms of the w-field. It is argued that for this expansion the proper- 
time regularization scheme has to be imposed on the operator level rather than on an 
expression in terms of the one-particle eigenenergies. The latter treatment is plagued by 
the inexact assumption that the Euclidean Dirac Hamiltonian and its Hermitian conjugate 
can be diagonalized simultaneously. It is then evident that approaches relying on counting 
powers of the w-field in the one-particle eigenenergies are inappropriate. 

Using the expansion of the action up to second order and employing a parametrical 
description of the soliton profiles the repulsive character of the u meson is confirmed. In 
the presence of the u meson the soliton mass is enhanced by a few hundred MeV. 
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1. Introduction and Motivation 

In the past few years the picture of baryons as being sohton solutions in effective 
meson theories has proven to be quite successful |1|]. As an example we refer to the 
appealing explanation of the almost vanishing matrix element of the axial singlet current 
for the proton, i.e. the so-called "spin-puzzle" 0. 

The investigation of the baryon structure in soliton models was started off by the 
consideration of the large Nc (number of color degrees of freedom) limit |P] of Quantum 
Chromo Dynamics (QCD), the theory generally accepted to describe processes being 
subject to the strong interacting. In this limit QCD was shown to be equivalent to an 
effective theory of weakly [0{1/Nc)) interacting mesons. Witten furthermore conjectured 
that baryons emerge as (topological) solitons of this meson theory Q. At small energies 
(or large distances) this effective theory should only contain the low-lying meson states, 
e.g. the pseudoscalar would-be Goldstone boson of chiral symmetry {TT,K,r]). A further 
requirement when modeling QCD was to maintain its symmetries, especially the chiral 
symmetry and its spontaneous breaking. These ideas were brought to practice by the 
rediscovery of the Skyrme model by Adkins, Nappi and Witten 0. The Skyrme 
model represents an effective theory of the pseudoscalar pions (and kaons in the three 
flavor generalization) and admits topological soliton solutions. In the spirit of Witten's 
conjecture the baryon number current, B^, was identified with the topological current 0]. 

Soon after this, vector mesons were added to the scenery in a way consistent with 



chiral symmetry [|, |^, [1^]. This has turned out to be especially pleasing since it brought 
together the fruitful concept of vector meson dominance (VMD) and the soliton picture 
of baryons [0. In particular the isoscalar-vector w-meson was shown to play a crucial 
role. It is directly connected to the baryon current via the stationary condition which 
generically reads 

{D^'' + mlg^'')u. = gB^. (1.1) 

Here D^'^ denotes a kinetic operator [e.g. d^d'^ — g^'^dpd^), which depends on the specific 
form of the model and g is a, coupling constant. For the time component, /i = 0, eqn 
( |1 . 1| ) is rather a constraint than an equation of motion. VMD then implies to identify 
the isoscalar part of the electromagnetic current as Wl 



ern^ 



jLo = -Y^' (1-2) 

with e being the electric charges. Transforming (|1.1|) to momentum space it is straight- 



forward to verify that the isoscalar radius of the nucleon is given by |TT 

rLo = ^1 + -^ (1-3) 

mi 

where r'^ = J d^rr'^Bo denotes the intrinsic baryonic radius. The additional contribution 
Q/mi, which is completely due to VMD, is crucial in order to reproduce the experimental 
value for rj^Q ^ (0.86fm)^. 

More recently efforts have been made to derive the effective meson theory from QCD 
motivated models of the quark flavor dynamics. In this context the model of Nambu and 



*As the integral over D^'^tOu vanishes J/^q ^^ properly normalized. 



Jona-Lasinio (NJL) [|T^ has acquired special attention |T^. The associated Lagrangian 



Nj-l 



j=0 



2^' '^2 



iV|-l 



2^2 E 1^(4^'^^)' + (4^57.9) j (1.4) 

is chirally symmetric. Here q denotes the quark spinors, rhP the current quark mass 
matrix and Nf the number of flavor degrees of freedom. The latter will be taken to be 
two from now on. As the coupling constants Gi carry dimension (mass)~^ the model is 
not renormalizable. This requires regularization and introduces an additional parameter, 
the cut-off A. For a sufficiently large coupling Gi the model indeed exhibits spontaneous 
breaking of chiral symmetry |]13|. This feature is reflected by a non-vanishing quark 
condensate (gg). 

The Lagrangian (|1.4|) can be rewritten in terms of composite meson fields by functional 



integral techniques |]T5|. This procedure is referred to as bosonization. In order to make 



these functional integrals well behaved it is mandatory to impose Feynman boundary 
conditions or, equivalently, to perform these integrals in Euclidean space. In Euclidean 
space the time-coordinate is analytically continued to Xq = —1x4^ = —ir. It also implies to 
continue the time components of the vector (V^ ~ ql^lq) and axial vector {A^ ~ qinl^q) 
fields 

Vo -^ iVi, Aq -^ lAi (1.5) 

and consider r, V4 and ^4 as Hermitian quantities. In Euclidean space the bosonized 
action finally is given as the sum of a purely mesonic term 

Ar. = j d'x (^-^tr(sts - m°(E + ^t) + (m°)2) + -^t^iY.V^ + A^A>^)^ (1.6) 

and a fermion determinant 

^F = Tr \og{ipE) = Tr log {i^ + f + 75^" {Pr^ + Pl^.^)) (1.7) 

i.e. ^NJL = -^F+^^m- Here P^^ = (l±75)/2 are the projectors on right- and left-handed 
quark fields, respectively. The complex field S describes the scalar and pseudoscalar 
meson fields, Sij = S°'t^J2 and Pij = P"r"/2: S = S* + iP. There is no problem with 
treating Am as it obviously is analytical in the meson fields. This, however, is not the 
case for Ap. As a matter of fact ^f is (in Euclidean space) not a real quantity and thus 
may be decomposed into real (Ar) and imaginary (Ai) parts 



Ar = -Tr log i^pEpE ' j , -^1=2^^ log [pEipE 

The above mentioned ultraviolet divergencies are completely contained in ArB Again 
the treatment in Euclidean space proves to be pertinent since the operator PePe'^ is 



^In order to obtain a finite value for Ai a suitable definition of the functional trace is mandatory. See 
section 3 for details. 



positive definite. Thus the proper-time regularization [|T^ prescription may be apphed 
to Ar resulting in the replacement 

l^ r'^ ds ( „ „ + 



Ar -^ --Tr y^ ^^ — exp \-spEpE ' j • (1-9) 

Expanding the fermion determinant in terms of the meson fields and their derivatives 



leads to an effective meson theory |15|. As this object is, by construction, a polynomial 
in the meson fields the analytic continuation back to Minkowski space is well defined. 
The resulting meson theory reasonably well describes the physics of the pseudoscalar and 
vector mesons. Furthermore this theory has many features in common with Skyrme type 
models which are supplemented by vector mesons. 

This effective theory also allows one to express physical quantities like the pion decay 
constant /vr in terms of the cut-off A and the constituent quark massB m = (E). In 
case the axial meson degrees of freedom are ignored this relation reads in the gradient 
expansion p 



/,= = i^^r 0,^ . (1.10) 




In practice the physical value f-^ = 93MeV will be employed to determine A for a given 
constituent quark mass m. 

Later on it was also demonstrated that the bosonized version of the NJL model con- 
tained static soliton solutions which extremize the Minkowski energy functional [1^. For 



the investigation of static field configurations the introduction of a Dirac Hamiltonian in 
Euclidean space, He, via 

ippE = -dr - hE (1.11) 

is useful because [dr,hE] = 0. The Euclidean energy functional is then expressed in 
terms of the eigenvalues of He- First calculations were constrained to the pseudoscalar 
mesons. Subsequently the isovector-vector [l^ and axial vector mesons |T^ were added. 



The associated field configurations have vanishing temporal components V4 = and 
A4 = rendering He an Hermitian operator. Thus the extraction of a Minkowski energy 
functional from the one in Euclidean space is straightforward. In this context it should be 
remarked that especially the inclusion of the axial vector meson provides further support 
for Skyrme type models because it turns out that then the valence quark orbit joins the 
Dirac sea and the baryon number is completely carried by the distorted Dirac vacuum 



Pi- 



The analytic structure of the energy functional changes drastically when the isoscalar 
vector tu-meson is included. In this case the soliton configuration involves a non-vanishing 
temporal component leading to a non-Hermitian operator He- So far three distinct ap- 
proaches to include the cj-meson into the static energy functional of the NJL model have 



been discussed in the literature |^-||2^. It is the main purpose of the present paper to 
illuminate the role of the tu-meson in the NJL soliton model and to study the behavior 
of the corresponding fermion determinant under transformations which mediate between 
Euclidean and Minkowski space. The above mentioned importance of the cj-meson for 
the description of static baryon properties requires a clarification of the situation. 



"''The difference m — ttiq — — 2Gi(gg) rcffccts tfie spontaneous brealcing of chiral symmetry. 



The remainder of this paper is organized as follows. In section 2 we investigate 
the analytical properties of the eigenvalues and -vectors of He- In section 3 we briefly 
review the derivation of the Euclidean energy functional and show that the analytic 
continuation which should transform this object to Minkowski space does not exist. In 
section 4 requirements on the Minkowski energy functional from a phenomenological 
point of view are derived. Furthermore several definitions for the Minkowski energy 
functional are discussed. Numerical simulations for the eigenvalues of He show that 
treating the cj-meson up to quadratic order represents a reasonable approximation. Using 
this result, it is argued that the same approximation for the Euclidean energy functional 
is pertinent. Then a Minkowski energy functional can actually be derived. The study 
of this approximation is presented in section 5. In that section we also demonstrate 
that counting powers of the cu-field in the eigenvalues of He leads to incorrect results for 
the module of the fermion determinant. As a matter fact, these misinterpretations are 

related to the appearance of Pe i^ ^^^ regularized version of the fermion determinant. 
Concluding remarks may be found in section 6. Technical details on the perturbative 
expansion used in section 5 are discussed in an appendix. 

2. Analytic Properties of the Dirac Hamiltonian 

In the introductory section we have argued that it is mandatory to analytically con- 
tinue forth and back from Minkowski to Euclidean spaces in order to obtain a meaningfully 
regularized energy functional. In the present section we will discuss the consequences of 
this continuation on the eigenvalues and -vectors of the static Dirac Hamiltonian He 
which can be extracted from eqn (|1.11|) . 



As inferred from eqn ( p.. 51 ) the time components of the (axial) vector meson play 
the key role for the study of the analytic properties. The ansdtze for the static meson 
field configurations are characterized by the fact that they commute with the grand spin 
operator G. This operator denotes the sum of total angular momentum (orbital plus 
spin) and isospin, i.e. G = I + cr/2 + t/2. Demanding in addition the proper behavior 
under parity the time component of the axial vector has to be zero while that of the vector 
field reduces to purely radial function uj{r). In Euclidean space the Dirac Hamiltonian 
may then be decomposed as 

hE = he + iuj{r). (2.1) 

Here h^ denotes the Hermitian part of He, i-e. h^ = ho,. Let us, for the time being, 
consider the simplest case when only pseudoscalar fields are present. The corresponding 
field configuration is assumed to be of hedgehog type 

S = m exp {ir ■ TQ{r)) . (2.2) 

Hence 

/i0 = a ■ p + m/3 (cos9 + i'^^.r ■ Tsin6) . (2.3) 

Later on we will also consider the case when p and ai fields are present as well. However, 
for the general discussion of the analytical properties the restriction to O and uj proves to 
be illuminating. For this discussion we generalize the Dirac Hamiltonian to an operator 
which depends on the complex variable z 

h{z) = /i0 + zuj. (2.4) 



The Euclidean space Hamiltonian, He then corresponds to z = i while the one in 
Minkowski space is associated with z = 1. We will refer to these values of z as Eu- 
clidean and Minkowski points, respectively. Except for the Minkowski point as well as 
z = —1, h{z) is non-Hermitian. In general we therefore have to distinguish between left 
{^u{z)) and right {^^{z)) eigenstates 

h{z)\m,{z)) = e,{z)\^!,{z)) (2.5) 

{i>M\h = e,{z){^>,{z)\ =^ h}{z)\^>,{z)) = e,{zY\^,{z)). 

Of course, the eigenvalues, e^^z), and -vectors, |\l'^(z)), parametrically depend on z. 
Adopting the set of free spherical grand spinors as basis |^ the matrix elements of h^ 

and uj turn out to be symmetric. Since h} {z) = h{z*) the eigenvalues and their complex 
conjugate are related via 

e^izY = e,iz*). (2.6) 

Furthermore a phase convention may by chosen such that the wave-functions are related 
by 

\^{z)) = |*(z)*) = |^(^*)). (2.7) 

For later discussions it will be useful to define real and imaginary parts of the eigenvalues 

e^i^, ^*) = I (e.(^) + e^z*)) eiiz, z*) = -'- {e^{z) - e^{z*)) (2.8) 

It is the main issue of this section to explore the 2;-dependence of the eigenvalues and 
-vectors once profile functions B(r) and uj{r) are made available. Of course, these test 
profiles can be chosen arbitrarily, however, they should (at least) be physically motivated. 
For convenience we will employ the profile functions displayed in figure 5.2 and constrain 
ourselves to the case m = 400MeV. The reader may consult section 5 on the origin of 
these profiles. Here it suffices to remark that the uj field satisfies an equation analogous to 
( |1 . 1| ) and is thus properly normalized. In order to study the behavior of the eigenvalues 
(and -vectors) in the complex plane defined by z we parametrize 

z = 6exp iiif) with < </? < 27r. (2.9) 

Obviously 5 = 1 and < (/) < 7r/2 parametrizes the continuation from Minkowski to 
Euclidean space and vice versa. However, we will consider 5 as a parameter in order to 
determine the radius of convergence in the complex 2;-plane. A priori it is not obvious 
that the energy eigenvalues exhibit any kind of analytic structure as these are roots of the 
characteristic polynomial (See e.g. appendix C of ref.|^.). This polynomial is in principle 



of infinite degree and only the restriction to a finite basis renders this degree finite as well. 
In figure 2.1 the dependence on 5 for two different energy eigenvalues is shown. On the 
left the energy of the valence quark orbit is displayed. This orbit is, by definition, the one 
with the smallest absolute value of the real part of the energy eigenvalue. I.e. the valence 
quark state is associated with the smallest |e^| (see also section 3). The valence quark 
state is distinct since in the presence of static meson fields it is the one which deviates 
most strongly from a free quark state. Obviously the valence quark energy represents a 
smooth function of (^ for 5 <1 while for 5 = 2 the derivative de^^i/dip becomes singular. 



At the "edges" y? assumes the values ±0.327r. These singularities are, however, not the 
consequence of a non-analytic behavior of the roots of the characteristic polynomial but 
rather a level crossing appears. For \ip\ < 0.327r a state with negligible yj-dependence 
carries the smallest |e^|. For \ip\ > 0.327r the role of the valence quark is then taken 
over by an orbit which strongly depends on if. Numerically we have confirmed that for 
the test profiles under consideration such level crossings in the valence quark channel are 
avoided as long as 5 < 1.3. This radius of convergence suggests that for the relevant 
parameter space {6 = 1) the eigenvalues of the Dirac Hamiltonian are analytic functions 
of z. This result was already obtained previously p6[. This is furthermore supported 



by considering other levels than the valence quark state. These actually possess a less 
pronounced dependence on u and have thus a larger radius of convergence. As an example 
the energy of the state with the smallest |e^| in the channel with grand spin and parity 
Qn _ 2+ ig displayed in figure 2.1. For this energy the radius of convergence is obviously 
larger than two. 

There is actually a much more elegant technique to investigate the analytical structure 
of the eigenvalues and especially to extract the radius of convergence. As these eigenvalues 
are supposedly analytic in z the corresponding Laurent expansion 

oo 

e.W= E c,,„(^o)(^-^or (2.10) 

n=— oo 

should exist. Here zq serves as the center of the expansion and may refer to any point in 
the complex plane. The Laurent coefficients c,^^n{zo) are defined by the Cauchy integrals 



2-ni JdA [Q - Zq) 



c.n{zo) = ^i, ^r''':i, (2.11) 



with Zq being in the interior of the integration contour dA, i.e. zq G A. Employing the 
parametrization ( = 5e*'^ + zq these Cauchy integrals become ordinary integrals 



CA^o) = 2^^, 



— / " e^ (de'^f + Zq) exp {-irup) dip. (2.12) 

ro" Jo ^ ' 



The computation of the integrals (2.12) obviously requires the eigenvalues of 



/i0 + (5e'^ + ^o) ^ (2.13) 

in the interval < (/) < 27r. Analyticity is then equivalent to the fact that the coefficients 
Ci/,n(-2^o) do actually not depend on b. On the other hand the Cy,ri(zo) are useful tools to 
determine the radius of convergence at Zq. For this purpose we start off with a small value 
for 5, numerically compute the eigenvalues of ( |2.13| ) and subsequently the integrals ( |2.12| ). 
Next the calculations are repeated assuming a somewhat larger value for 5. As long as 5 
is sufficiently small the Cu,n{zo) are indeed independent of 6. However, as a critical value 
is exceeded the Cu^nizo) vary with 6. This critical value has to be identified as the radius 
of convergence at the point Zq. For the valence quark orbit the results of this calculation 
are displayed in table 2.1. Again we have employed the test profiles which have been 
mentioned above. For the interpretation of this table (and similar once which will follow) 
it should be stressed that the entry "0.00" means that this value is zero up to the given 
accuracy. On the other hand the entry "0" implies that the corresponding coefficient 
vanishes identically up to numerical uncertainties which usually are of the order 10~^. 



Table 2.1: The coefficients Cvai,n(-2o) in the Laurent series (|2.10|) for the valence quark 
energy in units of the constituent quark mass. 
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Table 2.2: Same as table 2.1 for the state with the smallest (positive) real part of the 
energy eigenvalue in the G^ = 2"*" channel. All those coefficients which are not displayed 
either vanish [n < 0) or are negligibly small [n > 2). 



n 
n 



0.5 



^0 = 
1.0 



1.5 



1.1449 1.1449 1.1449 
0.0022 0.0022 0.0022 



0.5 



Zq = 1 

1.0 



1.5 



1.1453-1-0.00212 1.1453-^0.00212 1.1453-f0.002l2 
0.0020-0.0005Z 0.0020-0. 0005Z 0.0020-0.00052 



Evaluating the Laurent coefficients associated with an expansion around the origin of the 
complex plane we conclude that the radius of convergence for the energy of the valence 
quark energy eigenvalue lies in between 1.2 and 1.5. In the region of convergence the 
Laurent series actually reduces to a Taylor series without any poles. Then the analytic 
continuation from Euclidean to Minkowski space is unique for this eigenvalue. However, 
as will be discussed in the next section, the physical relevant situation is somewhat 
different. The fermion determinant is defined in Euclidean space only. Thus one has to 
perform a Laurent expansion around Zq = i. The results for the coefficients stemming 
from this calculation are also shown in table 2.1 and demonstrate that for this expansion 
the radius of convergence is less than 1.2. Taking into account that the distance between 
the Euclidean and Minkowski points is ^/2 we see that it is not possible to carry out the 
required continuation for the valence quark level, at least for the test proffies considered 
here. Again this non-analyticity originates from a level crossing as discussed above. Thus 
the identification of the valence quark level does not necessarily symbolize an analytic 
operation. For orbits other than the valence quark the continuation is less problematic 
because their response to the a;-profile is significantly weaker. As an example we have 
again considered the lowest state in the 2"*" channel. As can be inferred from table 2.2 
the dependence of the corresponding energy eigenvalue on z is almost negligible. 

There is one more conclusion which can be drawn from the results listed in tables 
2.1 and 2.2. In the region of convergence only those Laurent coefficients are sizable 
which correspond to a constant or linear dependence of the eigenvalues on z. In addition 
a small but non- vanishing quadratic piece may be present. This observation may be 



translated into the statement that the functional dependence of the energy eigenvalues 
on uj is at most quadratic. Of course, this again depends on the test profiles. However, 
our computations indicate that this is a generic feature for cu-meson profiles which are in 
some way related to the quark baryon current and hence are properly normalized. 

We have performed similar investigations not only for the eigenvalues but also for 
the eigenstates of h{z) (see eqn (|2.5|)). The symmetry of h^ and ui in the standard 
basis straightforwardly relates left and right eigenstates by complex conjugation. It is 
thus sufficient to only consider the behavior of the right eigenstates in the complex z 
plane. For the valence quark wave-function we display the dependence on z = 5ex]i{iip) 
in figure 2.2. Again we observe that the dependence is analytical for 5 < 1.0 and that 
for 5 as large as 2.0 the identification of the valence quark level exhibits an arbitrariness. 
Although the somewhat shaky curve for the lower component seems to suggest some kind 
of non-analyticity even for 5 = 1.0 this is actually not the case as we have verified by 
computing the Laurent series in analogy to ( p.lO|) . 

At the end of this section we thus arrive at the statement that except for the identifi- 
cation of the valence quark state the diagonalization of h{z) results in analytic expressions 
for its eigenvalues and -functions. For the identification of the valence quark state the 
radius of convergence is comparable to the distance between Euclidean and Minkowski 
points and depends on the center of the Laurent expansion. The identification will not 
cause problems if the self-consistent w-meson profiles turn out to be less pronounced than 
the test profiles considered here. Having computed a self-consistent solution it will be 
necessary to verify that the corresponding profile functions do not cause any level crossing 
for the valence quark state and the associated non-analytic behavior is avoided on the 
path from Euclidean to Minkowski spaces. Otherwise the solution has to discarded. 



3. Analytic "Properties" of the Energy Functional 

In the previous section we have examined the analytic properties of the eigenvalues and 
-states of the one particle Dirac Hamiltonian ( |2.4|) in the complex z-plane. Here we will 
investigate how (or whether at all) these properties propagate to the energy functional. 
For this purpose we will briefly review the derivation of the energy functional in Euclidean 
space and also emphasize on an important issue when evaluating the imaginary part of 
the action, Ai, see eqn ( pTSl) . 

Ignoring for the moment the regularization of the fermion determinant associated with 



the Hamiltonian (p.4|) it can formally be written as [27, 28 



Det {Pe) = CiT)exp{-TEo) J2 e^P -TNc J2 Vue. ■ (3.1) 

Here T denotes the Euclidean time interval under consideration. C(T) refers to a constant 
of proportionality which does not depend on the eigenvalues e^. Hence its explicit form 
is of no further relevance for our present investigations. The sum goes over all possible 
configurations {tiv\ of (anti) quark occupation numbers ri^ = 0, 1. The quantities ejy are 
defined in terms of the real and imaginary parts of the eigenvalues of e^, (see eqn (|2.8|) ) 

e, = |e^| + 2 sgn(e^)e^ = sgn(e^)e,. (3.2) 

The limit T ^ oo projects out the vacuum energy 

Eo = -^T.^. (3.3) 

from the exponential. Eqn. ( |3.1| ) provides a natural definition of the valence quark part 
of the energy 

Ev = E^ + lE'y = NcY.r^Ae'^\ + ziY^; ^r?,sgn(e^)e^. (3.4) 



It should be noted that the unregularized fermion determinant (|3.1| ) actually is a function 
of the energy eigenvalues e^ = e^ + iel rather than of their real and imaginary parts 
separately. Thus the complete fermion determinant [i.e. when the sum over all sets 
of occupation numbers is carried out) would represent an analytic function of z if it 
were finite and as long as the eigenvalues were analytic functions of z. For practical 
calculations, however, one considers only configurations with a definite baryon number 

B = Y.{vu-l)sgn{e^). (3.5) 

This corresponds to restrict oneself to a special set of occupation numbers without car- 
rying out the associated sum. Then an explicit dependence on a subset of the quantities 
e,y will induce a non-analyticity if sgn(e^(z)) = sgn{e^{z) + eiy{z*)) changes along the 
path connecting Euclidean and Minkowski spaces. From the discussion in the previous 
section it should be clear that such a change is most likely to happen for the valence 
quark orbit. All other states possess a sizable |e^| which is only slightly affected by the 
tu-meson. Let us for the moment consider the case with baryon number B = 1. The 
corresponding field configurations usually admit one distinct quark level which is referred 

10 



to as the valence quark orbit. The static energy corresponding to a configuration where 
this state is occupied is hence given by 

E = Eo + Ncv..ie..i (3.6) 

Nc ^ -. / r 



Y, ^u + Nc f r/val - - j Cval- (3.7) 



As already mentioned above, the states other than the valence state are almost unaffected 
by the w-field. In particular only sgn(e:J?^i) might vary when performing the analytic 
continuation. According to ( |3.5| ) the restriction to i? = 1 demands for the valence quark 
occupation number 

^vai = ^(l + sgn(e^,i)). (3.8) 

Substituting this result into eqn ( |3.7| ) yields 



E = ^\e^,,- Y: eJ (3.9) 

which turns out to be analytic in z in account on the above assertions on sgn(e^_^ya^i). 

The situation gets more involved when the fermion determinant is regularized. As 
already mentioned in the introduction the proper-time regularization requires to treat 
real and imaginary parts ( |1.8| ) separately. This leads to distinct treatments of ey{z) and 
ey{z)* . In turn the energy functional will depend on z and its complex conjugate, i.e.. 
E = E{z,z*). As only the vacuum part Eo undergoes regularization the transition from 
eqn ( |3.6| ) to eqn (|3.7| ) is, in general, prohibited. Thus there are (at least) two points which 
may cause the regularized energy functional to be non-analytic in z. Let us proceed step 
by step in deriving the regularized energy functional to see how this comes about. 

At this point we wish to apply the proper-time regularization scheme to the real part 
of the fermion determinant 

An=^ JTr log( -dr + h{z)) + Tr log(9. + h^z))^ (3.10) 

The functional trace is carried out by summing over the eigenvaluegil fi^ and e,y{z) corre- 
sponding to the operators dr and h{z), respectively. As we intent to extract Eq from the 
expression (^TTUj) we need to consider the limit T -^ oo. Then the sum over Qn becomes 
a spectral integral 

Nn^ f^ dv 



Z J-oo In ^ 
= f^r|E>og{l- + e„r}. (3.11) 



To this expression the proper-time prescription may be applied because the argument of 
the logarithm is positive definite 



An = -^T ^ / exp {-. t;2 + \e,{z)\' + 2vei{z)] } (3.12) 



*ri„ — {2n + l)n/T are the Matsubara frequencies. The fermionic character of the quarks reqmres the 
associated eigenfunctions to assume anti-periodic boundary conditions on the Euchdean time interval T. 
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Here it has to be stressed that this spectral integral cannot be obtained from the original 
definition of the proper-time regularization for the real part ( |1.9[ ). Substituting ( |1 . 1 1| ) 
and (0) yields for (^ 



A 



R 



-It. 

2 



~ ds 

1/A2 S 



exp 



-dl + h{z)h}{z) + lh{z) - h}{z)\ dr \ 



(3.13) 



because [dT-,h{z)] = 0. The spectral integral ( |3.12| ) could only be obtained from ( |3.13| ) 
under the assumption that h and h< can be diagonalized simultaneously. At the Euclidean 
point this assumption is violated at quadratic order in u^. Technically the difference arises 
from the fact that the usual rules for the logarithm have been used in order to derive 
( 3.11 ). However, the original definition 



|1.9|) of the proper-time regularization of Ar 
does not maintain these rules for finite A because the regularization curbs the Hilbert 
space. Stated otherwise: Different regularization prescriptions induce errors in the energy 
functional at order u;|. It should furthermore be stressed that in deriving eqn (|3.11 ) one 
has treated infinite quantities. This is avoided when the proper-time prescription is 



imposed at the operator level ( |1.9| , p.l3| ). In section 5 we will argue from the point of the 
numerical results that (p.l3|) represents the more appropriate starting point. Nevertheless 
let us continue with the investigation of the expression ( |3.12| ) in particular to later on 
compare the various approaches for treating the tu-meson in the NJL model. 

This spectral integral ( |3.12| ) converges absolutely permitting in particular to shift the 
integration variable v ^ v — el{z) 



Ar 



N, 



c 



T 



— / / — 6xp 

-oo 271 ^ Jl/h? S 



V 



M( 



(3.14) 



Performing the Gaussian integral allows one to extract {Ar -^ —TEq 



E^iz,z* 



Nc 



^E 



ds 



/A2 ^/A' 



ns'^ 



exp 



euiz) +eJz* 



(3.15) 



which obviously depends on z and its complex conjugate. In terms of the eigenvalues Qr. 
and ei,{z) the imaginary part of the action ( |1.8| ) reads 



^ n,u ^' 'n ^v\^ 



(3.16) 



Now it is important to note that the associated spectral integral 



IV - eAz) 



^ dv^ 
— log- 

oo 27r IV — ei^iz" 



(3.17) 



is not absolutely convergent. If it were one could easily eliminate the imaginary parts 
el{z,z*) of the energy eigenvalues from the integral (|3.17|) by reversing the sign of v in 
the denominator and subsequently shifting v by el{z, z*). We therefore need an extended 
definition of Aj. We will later on see that the principle value prescription 



2 ^ j-oo 27r IV — euiz* 



(3.18) 
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yields the desired physical symmetries such as the current field identity. The integral 
( |3.18| ) can now be performed using standard means 



Ai = -^Y.TV - dX 



V 



2 ^ J-oo 27r 7-1 i(v - Aef,) - e^ 



= ^^^T hm r'-^f'dX^ :'^^ ,. (3.19) 

Next the shift in the integration variable v — Xe^ — ^ t; is performed. It can be shown that 
the associated shift in the boundaries does not contribute as A^ — ^ oo. This yields 

^, = ^^ E /' ^A T hm r ^^. (3.20) 

Now the integral over the parameter A may be done. Although Ai is finite in the principle 
value formulation, the proper time regularization may be imposed by expressing the 
integrand as a parameter integral: 

which does obviously not diverge as A -^ cxd. We may finally extract the imaginary part 
of the Euclidean energy E^ from Ai —>■ —ITEq 

EU^, zl = ^ Ei^.iz) - 6.(^*))sgn(6.(^) + e^{z*)) 

J 1, Ai not regularized fioo^ 

\N',y{z,z*), Ai regularized ^ ' ' 

where we again made explicit the dependence on the complex variable z and its complex 
conjugate. Furthermore 



^,,..HJ.ri pii^ll^e* 



^Az) +eJz*) 



2A 



(3.23) 



are the vacuum "occupation numbers" in the proper time regularization scheme. The 
upper case in eqn ( |3.22| ) corresponds to the limit A ^ cxo. Previously p3| , |2^ the principle 



value description was not stated explicitly. Fortunately, due to the A integration the 
asymmetries in the boundaries cancel. Thus the final result remains unchanged. 

Analogously to the studies in section 2 we have then considered the total energy 
functional originating from the fermion determinant 

Edetiz, z*) = E^{z, z*) + iEl,{z, z*) + Ncr]..\t..x{z, z*) (3.24) 

in the complex z-plane. The imaginary part E^ is taken in its regularized form. The 
behavior of E{z, z*) corresponding to the parametrization ( p.9|) is displayed in figure 3.1. 
The cuts which occur for 5 = 2 stem from level crossings as described in the previous 
section. Otherwise E{z^ z*) exhibits a smooth behavior on the angle Lp. Although the 
valence and sea parts are treated differently the above indicated non-analyticity (see the 
discussion starting at eqn (|3.6|) ) due to the change of sgn f e:J^i(z, z*)\ is mitigated. This is 
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Table 3.1: Coefficients Cdct,n{zo) in the Laurent series for the contribution of the fermion 
determinant to the total energy (|3.26|) . These coefficients are measured in GeV. Also 
shown are the associated continuations to the Minkowski point Eu = Edct{z = 1) as 
obtained from the coefficients quoted. 







Zq ~- 


= 






zo 


= 2 




6 


0.1 


1.0 


1.2 


1.5 


0.1 


1.0 


1.2 


1.5 


n = -4 


0.000 


-0.007 


0.036 


-0.029 


0.000 


-0.008 


0.061-0.0852 


0.031+0.0432 


n = -3 


0.000 


-0.006 


-0.014 


0.135 


0.000 


-0. 004-^0.0322 


0.066-0.0132 


0.207+0.2452 


n = -2 


0.000 


-0.064 


-0.136 


-0.260 


0.000 


-0.045+0.0122 


-0.051-0.0982 


-0.111-0.1922 


n = -1 


0.000 


-0.005 


0.069 


0.330 


0.003Z 


0.015+0.312? 


0.039+0.3702 


0.058+0.636Z 


n = 


1.208 


1.348 


1.425 


1.572 


0.706-f0.6552 


0.922+0.6542 


1.024+0.5522 


1.202+0.4812 


n = 1 


0.650 


0.655 


0.656 


0.618 


0.689-^1.1642 


0.675+1.1552 


0.649+1.0782 


0.602+1.0692 


n = 2 


0.740 


0.742 


0.744 


0.717 


0.484-0. 024z 


0.528-0.0132 


0.512-0.0682 


0.510-0.0482 


n = 3 


-0.014 


-0.013 


-0.012 


-0.032 


0.008-f 0.1642 


0.004+0.1242 


-0.028+0.0882 


-0.033+0.1112 


n = A 


0.042 


0.041 


0.039 


0.028 


0.075-0.0542 


0.023-0.0032 


-0.003-0.0202 


0.009+0.0042 


EuiGeV) 


2.627 


2.700 


2.673 


2.782 


2.5234-0.0352 


2.712-0.0362 


2.711+0.1372 


2.715+0.0292 



a special feature of the proper-time regularization which contains the equivalence between 
the limits e^ — > and A — > cxd in eqns (|3.15|J3.2^ ). This is intuitively clear because these 
regulator functions only depend on the dimensionless ratios e^/A. The smooth behavior 
of E{z, z*) in the complex 2-plane for the physically relevant radius S ^ 1 suggests that 
E{z, z*) resembles an analytic function. This, of course, would be somewhat astonishing 
because z* cannot be expressed as a power series in z without trading in an essential 
singularity. As E{z,z*) depends on both z and z* this essential singularity should carry 
through to the energy functional. In order to clarify the analytic properties of the energy 
functional (|3.24|) we have also examined the corresponding Laurent series 



Edet{z) ■= Yl Cdet,n(^o)(^ 



^oj 



(3.25) 



Again we considered expansion centers located at the origin of the complex plane {zq = 0) 
and the Euclidean point [zq = i). The results are summarized in table 3.1. 

There are several observations to be made. First, we find non-vanishing coefficients 
Cdot.n with n < even for parametrizations which led to non-singular one-particle energies 
eu{z). Second, these coefficients obviously depend on the path (dA) along which the 
Cauchy integrals 



Cdet,n(^o) = TT- f d(—^ 
27r2 JdA (C — 



(C,C 



(C - ^o) 



n+1 



(3.26) 



are computed. The path dA is again parametrized by ^ = 6exp{iip) with (0 < ^9 < 27r). 
Even for small areas A which include the expansion center {zq G A) the coefficients Cdet.n 
do not converge. Especially the coefficient of the constant term Cdet,o varies by about 
20% when decreasing 6 from unity to 0.1. For both values of 6 the one-particle energies 
were found to be analytic, see section 2. Furthermore we extract from our numerical 
studies that (for a given path dA) the coefficients Cdet,n decrease only very slowly as 
\n\ gets larger. Thus we have accumulated (at least) three arguments which support 
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the conclusion that the Laurent series ( p.25| ) does not converge: The coefficients Cdet.n 
with n < do not vanish, the radius of convergence shrinks to zero, and the Minkowski 
extracted by analytic continuation has a non- vanishing imaginary part. We have then 
considered the continuation to the Minkowski point Em = E^ctiz = 1). Again we find a 
dependence on 6. This dependence is somewhat weaker than in the case of the coefficients 
Cdet,n- We should, however, note that for the computation of Em only those coefficients 
which are displayed in table 3.1 were taken into consideration; i.e. \n\ < 4. The reason 
is that the numerical computation of the coefficients with larger n contains significant 
errors for small M. A further evidence for the non-existence of the series (|3.25|) can be 
deduced from the fact that Em is not real when the expansion is centered at the Euclidean 
point. We close our considerations on the structure of E^et with the statement that the 
associated Laurent series has vanishing radius of convergence. Stated more drastically, 
the analytic continuation of E^ct does not exist. And in general 

E^,,{z,z*)^E^et{z). (3.27) 

The Euclidean point resembles the only exception to this inequality. 

The non-existence of an analytic continuation for 5 < 1 is purely caused by the 
proper-time regularization of the Euclidean action. As this regularization treats real and 
imaginary parts of the one-particle energies differently it is interesting to consider other 
regularization schemes which do not have this property. The most suggestive approach 
would be to regularize the vacuum energy ( p.3|) by employing a sharp cut-off and include 
all states with |e^| < A in the summation. We have also studied the analytic properties 
of an energy functional defined in this way. Indeed it yields a small but finite radius of 
convergence. In this area of convergence the coefficients Cdet.n of negative index n < 
vanish. Convergence is lost as one of the real parts of the energy eigenvalues crosses ±A 
along the path connecting Euclidean and Minkowski space. Since numerically this radius 
of convergence is found to be significantly smaller than unity (which is the physically 
relevant one) we will not further pursue this treatment of the energy functional. Further- 
more a sharp cut-off regularization should be discarded because it does not even preserve 
global gauge invariance. 

4- Consistency Conditions for the Energy Functional 

In the last section we have learned that no analytic continuation of the regularized 
energy functional exists. One therefore has to reside on physical arguments, as e.g. gauge 
symmetry, in order to define (or motivate) an energy functional in Minkowski space. This 
will be the contents of the first part of this section. Actually three distinct approaches 
to the Minkowski energy functional are known in the literature. Historically the first 
to consider the a; meson in the context of the NJL model were Watabe and Toki 



Secondly a different calculation performed by researchers at Bochum University ^^ ^ 



while the third approach was put forward by the present authors |2^, |2^. In the second 
part of the present section we will examine these approaches from two points of view: 
(1) their relation to the Euclidean energy functional and (2) their consistency with the 
conditions which will be derived in the first part of this section. 

A constant shift of the tu-field a;(x) -^ uj{x) + ujc is equivalent to a global gauge 
transformation and should leave the fermion determinant unchanged, at least for topo- 
logically trivial configurations. The non-regularized form of the fermion determinant 



^The Cdot,n's are of the generic form / x (5 ", c/. eqn ( ^.12| ). Thus numerical errors in the integral / 
are multiplied for small 5 and large n. 
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np can easily be seen to possess this property. This is mostly a consequence of the 
principle value prescription for the imaginary part of the action ( |3.18| ). At the Euclidean 
point a constant shift in u) causes the imaginary part of the energy eigenvalues to change 
accordingly 

el^^el + ujc as uj{x) —^ uj{x) + ujc- (4.1) 

This straightforwardly transfers to that part of the energy which stems from the fermion 
determinant 

Eo^Eo-^^^cEsgn(ef) and Ey ^ Ey + iNcuo,Y.^,sgn{e^) (4.2) 

such that 

^det ^ ^dct + iNcBuj,. (4.3) 

This is just the desired result that shifting the lo field by a constant changes the energy 
(associated with the fermion determinant) by an amount proportional to the baryon 
number. This is not surprising since a constant lOc represents a chemical potential for 
the quarks. Since the unregularized form of the fermion determinant has been seen to 
be analytic it is obvious that the relation ( |4.3|) can be continued to Minkowski space. 
I.e. the Minkowski energy changes by NcBuj,.. We therefore require a suitable definition 
of the regularized energy functional in Minkowski space to transform accordingly under 
u}{x) — > ui{x) + ijJc- At this point we also recognize the reason why a sharp cut-off 
regularization in Minkowski space does not meet this requirement. In case lOc shifts one 
of the eigenvalues across this sharp cut-off this special state is dropped in the sum (|3.3|) . 
Then the transformation property (O) is no longer valid. The same behavior can be 



observed for other regularization schemes, as e.g. a three dimensional momentum cut- 
off, which operate in Minkowski space. If one wanted to implement the transformation 



relation (O) for a three dimensional momentum cut-off in Minkowski space one would 
inevitably be led to a regularization function which depended on lOc- This feature is, 
of course, undesired. In Euclidean space, however, only e^ enters the regularization 
condition. Thus the sharp cut-off or three momentum regularization in Euclidean space 
is consistent with ([4.3| ) because the shift ( ^?1| ) only effects the imaginary parts e^. 

It should be remarked that the transformation property ( [4.3| ) holds in Skyrme type 
models as well. In those models the coupling between to and the chiral field is of the form 



NcUJ^B^ |]5^. Shifting the tu field by a constant leads to the variation of the energytl 

Nc I dh iu.Bo = Ncto, I dh Bo = Ncto.B. (4.4) 

In order to extract the baryon current one introduces an external space-dependent 
gauge field bf^{x) for the f/y(l) symmetry. The baryon current B^ is then identified as 
the coupling linear in bf^{x), i.e. 



., _ 5A[b] 



(4.5) 



* Other couplings appear as total derivatives of topologically trival objects and thus do not contribute 
when shifting the w-field. 
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The action A[b] is defined by replacing y^with Y+ ^ in gQh ( p. .71 ). This guarantees the 
local f/y(l) symmetry. Under the functional integral the vector meson fields may be 
shifted by an arbitrary amount (again up to topological defects). This allows one to 
transfer the dependence on the gauge field to the mesonic part of the action ( |1.6| ) and 
straightforwardly obtain the baryon current 



i.e. the current field identity ( |1.2| ) holds. In chiral soliton models the identity 
commonly is a consequence of the stationary condition for the cj-field. This condition 
is obtained by extremizing the Minkowski energy functional. Thus we have at hand a 
second consistency condition for the Minkowski energy functional. For the static a;-field 
this can also be interpreted as the normalization of the cu-field 

I d^ruj{r) = G2B. (4.7) 

We will require the equation of motion for uj{r) to provide this normalization. This 
normalization also has the important consequence that the inclusion of the u meson into 
the soliton calculation does not necessarily lower the static energy although the number 
of degrees of freedom is increased. 

The first approach to incorporate the uj meson in the NJL soliton calculation was 
carried out by Watabe and Toki|2^. Their calculation was completely performed in 
Minkowski space. In this treatment the contribution of the fermion determinant to the 
energy functional has been defined as 

4el = E^o' + E^' = ^T.r., -^ -P h4) + iVcr^vaievai (4.8) 

z ^ J1/A2 v47rs'^ ^ ^ 



wherein e^, denote the eigenvalues of h(z = 1) (|2.4|) . Although this is formally the analytic 
continuation of ( p.l5| ) it has the disadvantage that the above stated consistency conditions 
are not satisfied. E.g. the shift uj{x) ^ uj{x) + Uc leads to e^ ^* e^ + uJc and thus all 
powers of lOc appear in ( [4.8| ) after the shift. HO|Wever, the requirement ([4.3|) only allows 
for a linear term. Furthermore as the derivativeP 



dE^^^ Nc 



de^ 



E ^r^ = -^ E sgn (6,) erfc -f U -^ E sgn (e,) (4.9) 



A 



N, 



2 . 



the normalization condition (|4.7| ) is not satisfied. Since the analytic continuation of (|3.24| ) 
does actually not exist it is not surprising that the formal continuation does not meet the 
consistency conditions. As these are not met by (|4.8| ) we do not consider it as a suitable 
definition for the Minkowski energy functional. 



Soon after this, Schiiren et al. |21, 22|| presented an approach which assumes as 



starting point the Euclidean action in the proper-time regularization scheme ( |3.15|J3.2^) . 
In order to reach the Minkowski point these authors "interpreted" the eigenvalues e^ of 
the operators h{z = ±1) as the continuations of ej,(±i). Then 



^i2 = f E {£ ^exp (-!(.: + .,f) H- i(.: - .;)sg„(.: -H .;)} (4.10) 



^This derivative can, however, be identified as the regularized baryon number carried by the polarized 
Dirac sea. 
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was defined as the contribution of tlie Dirac sea to tlie energy functional^. Taking ana- 
lyticity of tlie one-particle energy eigenvalues for granted eqn ( [4.10| ) amounts to 



.(2) 



E^:( = Edet(l,-l). (4.11) 

Apparently this interpretation cannot be characterized as an analytical continuation. 
Advantageously the definition ( [4.10| ) allows one to identify certain powers of the tu-field 



in the sence that e^ is even while e^, is odd in agreement with e^ and e^, respectively. 
However, already the quadratic order in u (contained in e^) of E\J^ differs significantly 

from the expansion in uj^ which starts from the expression PePe- This subject will 
be discussed at length in section 5. Here it is merely important to keep in mind that 
counting powers of oj in e^ is (at least) suspicious; most likely faulty. In contrast to 
the work of Watabe and Toki [^ the treatment of Schiiren et al. satisfies the above 
derived consistency conditions. Shifting the u;-field by a constant implies e~ ^ e~ + uOc 
while e^ remains unchanged. Substituting this transformation prescription into eqn (|4.1(]| ) 



immediately provides the analytical continuation of (|4.3| ). The normalization condition 
(^4.7| ) is verified by integrating eqn (4.37) of ref. p^, which represents the equation of 



motion in this approach, over the whole coordinate space. 

The approach of the present authors to include the uj degree of freedom into the NJL 
soliton calculations is based on an observation made in section 2. When investigating 
the analytic structure of the single quark eigenenergies it was found that in the physical 
relevant region the Laurent series centered at the Euclidean point contained only two 
significant terms, a constant and the expression linear in the complex variable z, see 
tables 2.1 and 2.2. Transferring this approximation to the total fermion determinant 
implies to identify 

Edet(^) ~ E^{i, -i) + zE'^{i, -i) (4.12) 

because at the Euclidean point the real (imaginary) part contains only even (odd) terms 
in uo^. This approximation is furthermore supported by the fact that the ambiguities for 
regularizing Ar are of order ul (see the discussion proceeding eqn ( p.l3| )). By construc- 
tion, the approximation ( ^.12| ) represents an analytic function of z which can be continued 
to Minkowski space yielding 

E'^2t = E^{h-i)+Ei{i,-i). (4.13) 

(3) 

It should be stressed that in contrast to the two other approaches the evaluation of £'jg( 
indeed requires the diagonalization of a non-Hermitian Dirac Hamiltonian. The one- 
particle energies which are evaluated in Euclidean space enter this expression for the 
Minkowski energy functional whence it is clear that the condition ([4.3|) stemming from 
the global gauge invariance is satisfied. The normalization condition (^4.7|) is proven to 
hold by integrating the equation of motion for the uj profile function (Eqn (B.6) of ref . ||24|| ) 
over the coordinate space!. 

To conclude this section we would like to make some remarks on a possible saddle point 
approximation to the functional integral in Euclidean space. In this approximation the 
stationary conditions for the profile functions are obtained by identifying the configuration 



■i-The imaginary part (3.22) has been taken in the unregularized version. 

§In ref. |24| we have considered the imaginary part of the action in the regularized as weU as in the 
unregularized form. The consistency conditions are only satisfied for the unregularized form. 
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with the smallest real part of the Euclidean action. The reason is that this configuration 
is assumed to dominate the functional integral. The resulting equation of motion acquires 
the form 

= ^-^§^ (4.14) 

Oip{rj 

where -Eij(i, —i) = E^{ii ~i) + Em{i) also includes a part which originates from the 
mesonic part of the action ( |1.6|) continued to Euclidean space. Furthermore ip denotes 
the meson field under consideration. As Eji{i, —i) contains only even powers of a; the 
saddle point equation (|4.14|) possesses the undesired solution u = 0. This solution does 



clearly not satisfy the global gauge symmetry condition ( [1.3| ). Thus the saddle point 
condition has to be amended in order to comply with the consistency conditions. Obvi- 
ously the extremization of the functional in which the fermion determinant is replaced 
by the approximation ( [4.12| ) at z = 1 represents a suitable extension of the saddle point 
equation. Although this replacement is not unique it provides further justification for our 
treatment of the u meson. 

5. Time Components up to Quadratic Order 

Up to this point the situation looks somewhat discouraging because we have merely 
accumulated arguments why several approaches cannot be considered correct. These 
approaches started on a complete non-perturbative treatment of the u meson in the NJL 
soliton calculation. We have seen that then an exact analytic continuation from Euclidean 
to Minkowski space is impossible for the regularized module of the fermion determinant. 
In this section we will thus return to a perturbative type of treatment for the cu-field and 
investigate up to which extent one can make contact to the non-perturbative calculations 
discussed in the previous sections. 

We start off with the Euclidean Dirac Hamiltonian h^ + iu^ and consider lo^ as a 
perturbation to h^. Standard perturbation techniquesci may be employed for the energy 
eigenvalues to illustrate the a;4 dependence 

e, = el + el + el... (5.1) 

with 

e^ = t{u\u4W) and e^ = - 2^ ^^^ • (5.2) 

Here e° denote the eigenvalues of Hq and \p) the corresponding eigenstates. The distinc- 
tion between real and imaginary parts is trivial 

R V- {T^\uJi\lj){li\uJi\u) / ■/ I I \ (K'i\ 

^u=^u-2^ ^Zr~o ^^^ ^v = ^\^\^4\^) ■ (5-3) 

In table 5.1 we summarize the numerical results of this perturbative calculation for the 
low-lying states in the G^ = 0"^ channel which also contains the valence quark level. 
Again we have adopted the standard reference profile functions already used to gain the 



* These techniques assume a complete Hilbert space. However, regularization effectively reduces this 
space. 
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Table 5.1: Perturbation for the one-particle energy eigenvalues of the complex Dirac 
Hamiltonian in the C^ = 0~^ channel. All numbers are in units of the constituent quark 
mass m. 





Perturbation 


Exact 


6° 


4 


4 


total 


hiz = z) 


-1.522 


0.1092 


-0.016 


-1.537-h0.109z 


-1.550+0. 119z 


-1.293 


0.093i 


-0.026 


-1.319+0m3i 


-1.333+0.0862 


-1.118 


0.069i 


-0.042 


-1.159+0.0672 


-1.150+0.0422 


-1.029 


0.018i 


-0.022 


-1.051+0.0182 


-1.037+0.0692 


0.382 


0.595i 


0.019 


0.401+0.5952 


0.397+0.6052 


1.071 


O.Olli 


0.006 


1.072+0.0112 


1.072+0.0092 


1.246 


0.034i 


0.004 


1.240+0.0342 


1.250+0.0292 


1.481 


0.056i 


0.009 


1.491+0.056Z 


1.490+0.0502 



results presented in sections 2 and 3. We observe that the perturbation series seems to 
converge quickly and that the second order (e^) is already negligibly small. Furthermore 
the comparison with the exact diagonalization of the Euclidean Dirac Hamiltonian shows 
reasonable agreement at the level of a few percent. We have observed a similar behavior 
for states not listed in table 5.1. This corroborates our earlier assumptions which led to 
the approximation ( |4.12| ). 

It is then intuitive to transfer this expansion to the complete energy functional at the 
Euclidean point. This will then be an analytic function in 004 which can straighforwardly 
be continued to Minkowski space via u^ — > —iuQ. Up to quadratic order in 004^ the energy 
functional in Euclidean space is given by 

TV / f° \ 

£^"^''[0,^4 = Eo[e] +2iVcr/vai(vaV4|val) -z^5:sgn(e°)erfc ( f j (ulu^W) (5.4) 






|(z/|ci;4|val)| A^( 



pO _ ,0 

val u 






m^ 



drr'^ul- 



Here ii^ol©] refers to the energy functional associated with the Hermitian part of the 
Dirac Hamiltonian Hq. Also the subtraction for the trivial vacuum configuration = 
is contained in i?o[0]- For the valence quark contributions (proportional to ?7vai) the 



application of ( |5.1| ) and ( |5.2| ) is straightforward. The corresponding calculation for the 
Dirac sea, which yields the regulator function 



I i sgn(^.)erfc(|^|)-sgn(£^)crfc(|^|) 

/ [eu, ef,; A) = <( 2 e^-e^ ' e,. ^ e^ 

0, e^ = e^ 



(5.5) 



may be found in the appendix. Here it is important to remark that this regulator function 
has vanishing diagonal elements (/ (e,y, e^; A) = 0) as a consequence of the second order 
perturbation formula (|5.2|) . However, the limit e^ -^ e^ is not assumed smoothly. The 
last term in eqn (|5.4| ) stems from (|1.6|) with the gradient expansion result [|15] G2 = 
3m^/2m^/^ substitutedQ. In case one wants to consider the model with the imaginary 



^^Here we do not include the effect of tt — ai mixing. 
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Table 5.2: Perturbation for the energy corresponding to the fermion determinant at the 
Euchdean point. All numbers are in MeV. 



Ev 

Eq + Ey 



Perturbation 
nth 1 St Olid 



458 714z 
1261 681z 



19 
-553 



pEucl 
-^2 



477+714Z 
724+681Z 



Exact 

Eiz = i,z* = —i) 



476+726Z 
705+692Z 



part of the action not being regularized the complementary error function in eqn ( |5.4| ) 
has to be replaced by unity. In that case the consistency condition ( |4.3|) is satisfied. 
This is easy to see since for a constant cj-field (z/|u;c|/i) = ojc^uti- As already mentioned 
the diagonal elements do not contribute at quadratic and higher ordeifl hence the whole 
change under the constant shift is proportional to the baryon number. From the numerical 
results listed in table 5.2 we find that the expansion ( |5.4| ) is reliable at the 3 percent level 
although the contribution from the Dirac sea quadratic in uo^ is unexpectedly sizable. 

Although these results are encouraging we will now argue that /(e;y,e^;A) in eqn 
( |5.5| ) is not the correct regularization function at quadratic order. We have already 
remarked the unpleasant feature that / (ei,,e^; A) is discontinuous as e^ — *• e^. A second 
disadvantage of / relates to the Minkowski energy and the resulting equation of motion 
for ujq 



2 / rfrV2/C(r,r>o(r') 



^val 




^-^ 1 / dml,{r)^>^^,{T). (5.6) 



Here we have defined the bilocal kernel 

47r mlfl 5{r - r') 
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dQ'^^li 



r')^,{r')f 



^°,eJ;A 



.t, 



M/l(r)M/, 



(5.7) 



Since \l'iy(r) = (r|z/) denote the coordinate representations of the ejgenstates |z/) of Hq 
this kernel functionally depends on the chiral angle G. Obviously /C(r, r') is symmetric 
in its arguments. The RHS of the equation of motion ( p.6| ) is, of course, nothing but 
the quark baryon density in the regularized formi. Upon integration of this equation one 
easily verifies the normalization condition ( |4.7] ). Here the important ingredient is again 
the fact that the regulator functions / has vanishing diagonal elements. By the virial 
theorem one can then eliminate the linear term in the Minkowski energy functional for 
profiles which satisfy (p.6|) 



E^'''''[e,u;o] 



En 



+ - drr / dr'r' )C{r,r') u!o{r)u!o{r') 



{51 



■l-This is due to the fact that the perturbation of the wave-function is orthogonal to the unperturbed 
state in Rayley-Schrodinger perturbation theory. 

§For our test profiles the question of regularizing the baryon number or not plays a minor role since 
numerically the regularized baryon number comes out to be 0.96. 
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Only in case that all eigenvalues of this kernel are positive the Minkowski energy is positive 
definite. For the numerical treatment we have discretized the radial coordinate leaving 
this kernel as a matrix ]C{r, r') -^ ICij. For various test profiles B we have found that this 
matrix possesses negative eigenvalues. This can easily be understood by considering / in 



eqn ( ^.5|) . The dominant contribution comes from states with e^ ~ e^. Assuming e^ > e^ 
(without loss of generality since / is symmetric) we see that the denominator is positive. 
On the other hand the numerator is negative because the complementary error function 
decreases monotonically with the argument. Thus / may assume large negative values. 
Stated differently, the fact that the Minkowski energy functional is not positive definite 
for profiles satisfying the stationary condition ( ^.6|) is a consequence of the discontinuity 
oi f {eu,ef,;A) as e^ -^ e^. 

We have thus seen that counting powers of the cj-field in the eigenvalues of the 
Hamiltonian h^ + iu^ leads to an ill-defined expression for the real part of the fermion 
determinant in the proper-time regularization ( |1.9|) . We also conjecture that the non- 
positiveness of the energy functional based on this power counting represents the reason 



why the authors of ref. ||2T[ were unable to find self-consistent solutions for physically 
relevant parameters. 

To resolve this problem let us next examine Ar in more detail. In order to expand 
Ar in terms of 004^ one starts off at the module 

pEpl = -dl + /i| + liuj^dr + % [/le, ^4] + ^l- (5.9) 

Performing subsequently an expansion of 



Ar = -TT\og[-d^ + h'Q + 2iuj4dr + i[he,uj4+ujij (5.10) 

1 roc (j^g , , 

-^ --Tr / — exp(-92 + /i| + 2zcj45^ + z[/ie,cj4]+cjn (5.11) 

corresponds to imposing the proper-time regularization scheme at the operator level. 
As has already been explained in section 3 the regularized form ( |3.15|) follows from an 



intermediate definition of the proper-time regularization and involves the assumption of 
simultaneous diagonalizability of h and h} . This is exactly the point where ci;4-power 
counting of the eigenvalues fails to provide the correct result for the expansion of Ar. 
We therefore expect different results for the action functional when adopting ( |5.11|) as 



starting point for expanding in uj/^ rather than expanding the energy eigenvalues. Note 
that the expression ( p.9|) is still exact in a;4, i.e. no expansion has been performed up 
to here. The vacuum contribution to the real part of the fermion determinant is then 
expanded up to quadratic order in uj^ by methods which have previously been worked out 
in the context of quantizing the chiral soliton [^ and has subsequently been extended 
to study fluctuations off the chiral soliton ||29[. We refer the interested reader to the 



appendix where the calculation is performed. The net result of this calculation is that 
the regulator function / deflned in eqn (|5.5|) has to substituted by 



., lsgn(€jerfc(^)-sgri(£„)erfc(i) A e-<<.™" - e-l<'/«' _„ 

^ ''- ^- ^' = 5 7-^^ ^ irrji ■ (5-12) 

The corresponding replacements have also to be performed in the equations (|5.4| ) and 
The latter replacement deflnes the bilocal kernel /C which then enters the new 
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equation of motion and energy functionals obtained by the corresponding substitutions 
in equations (^]6|) and (|5.8|), respectively. Apparently these two regulator functions (/ 
and /) only agree in the infinite cut-off limit, A ^ oo. This is not astonishing because 
in this limit the introduction of the Hermitian conjugate of the Dirac operator is not 
mandatory in order to (formally) evaluate the fermion determinant. It is also noteworthy 
that the regulator function / ( |5.12| ) is identical to the one which appears in the moment of 
inertia when the semi-classical cranking approach is applied to quantize the chiral soliton 



27| . This is intuitively clear because for that problem the perturbation (proportional to 
the isospin generators) in the Dirac Hamiltonian is also static and anti-Hermitian. In 
that context it has been shown that the limit 

lim / (e^, e,; A) = (5.13) 

is assumed continuously (see chapter 5 in ref. |^). Numerically we have also verified 
that the associated kernel /C possesses positive eigenvalues only. Thus the corresponding 
energy functional is positive definite for uj profiles which statisfy the modified equation 
of motion (/C inserted in eqn (|5.6|) ). Due to the limit (^.13|) it is clear that the Minkowski 



energy functional with / replaced by / satisfies the consistency conditions derived in 
section 4. For our standard test profile we find that the second order in ujq contributes 
about -129MeV to the Minkowski energy. This is significantly smaller (in magnitude) than 
for the regularization function / (see table 5.2). Thus the numerically results put further 
doubts on the validity of counting uj^ powers in the energy eigenvalues. On the other hand, 
employing the correct regulator function / not only satisfies our consistency conditions but 
also provides a quickly converging series. The zeroth and first order contributions may be 
read off from table 5.2. Thus we have, for the first time, available a well suited Minkowski 
energy functional for the NJL model which includes the u meson. For completeness and 
later reference we list it here again although the (formal) deviation from ( p.4|) is only 
minor 

E^'^'^lQ.Uo] = Eo[e]+iVcr/vai(vaVo|val)-^5:sgn(eO)erfcf -^ j(z/|^o|^^) (5.14) 
^ |(;/|^o|val)|^ A^c ^ / A^ 1/ I i \|2 ^vrmj/^ r ^ ^ 

For UJ profile which statisfy the associated stationary condition this can be decomposed 

as 

E^'^[Q,uJo] = Eo[e]+E4Q,uj] (5.15) 

where 

E^[<d^uj] = 77 / drr"^ / dr'r'^JC {r,r') uJo{r)uJo{r') (5.16) 

denotes the contribution due to the u field. 

In the next step one has to construct the self-consistent solution which extremizes 
the functional ( ^.14|) . Unfortunately this is not an easy task. The reason is that the 
kernel /C functionally depends on the chiral angle B. Thus the stationary condition 
for G requires the functional derivative of the wave-functions "^^Ir) with respect to G. 
Such a calculation goes beyond the scope of the present paper and we will pursue a 
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Table 5.3: Parameters describing the energy minimum as well the corresponding total 
energy in Minkowski space, see eqn ( p.l4| ). Also shown is the contribution of the uj field 



(|5l6|) 



m (MeV) 


^mm 


^Mink (^^gy) 


E^ (MeV) 


i?-- (fm) 


^Mmk (^gy) 


E^ (MeV) 


400 


0.852 


1586 


321 


0.713 


1590 


317 


450 


0.776 


1685 


380 


0.775 


1685 


378 


500 


0.729 


1769 


428 


0.830 


1768 


419 



different approach to find the energy minimum. First, we introduce parameters {x} to 
describe the shape of the profile function B(r) and evaluate the corresponding eigenvalues 
and -energies of h^. These are then functions of the parameters {x\- Iii a second 
step this allows to construct the kernel /C as a function of these parameters as well. 
Via the stationary condition for cuq, eqn ( p. 61 ) with / replaced by /, also ojq becomes a 
function of these parameters. As already mentioned in the introduction this stationary 
condition commonly represents a constraint rather than an equation of motion. This 
fact furthermore justifies our procedure to parametrize the chiral angle while solving 
exactly for the uj profile. Technically this stationary condition is an inhomogeneous 
integral equation which is solved by inverting the discretized bilocal kernel K,ij. Finally 
the parameters are tuned such that E"^™"^ acq^jirgg ^ minimum. Since the main effect 
of the UJ meson on the chiral angle is of repulsive character it is intuitive to introduce a 
parameter which allows to vary the spatial extension of B(r). In a first calculation we 
therefore introduce the breathing mode parameter x = -^ via 



eA(r) = es.c.(Ar) 



(5.17) 



where Os.c.(^) denotes the self-consistent soliton solution to the problem without the uj 
meson. The resulting energy is displayed on left of figure 5.1 as a function of A. Secondly 
we have adopted the description |^ 



e 



Rrr 




rriTrRr, 



exp{m^{Rra -r)) 



r <R„ 
r>R„ 



(5.18) 



for the chiral angle. Since Rr is chosen such that the derivative of G is continuous 
the chiral angle only depends on R^ parametrically. The results for the energy cor- 
responding from the parametrization ( |5.18| ) may also be found in figure 5.1. Actually 
these two parametrizations yield almost identical results. For a constituent quark mass 
m < 400MeV no local minimum is obtained and the trivial minimum {E = 3m) corre- 
sponding to three free quarks is assumed. In the latter case the chiral angle gets very 
narrow while the cu-field tends to zero except within a vicinity of the origin such that 
I J d^r uj\ > and | J d^r uj^\ ^ 0. As the constituent quark mass is increased the local 
minimum gets more pronounced. In table 5.3 we list the parameters describing the min- 
imum as well as the resulting total energy. There A™° and i?™™ denote the parameters 
associated with the energy minimum. As already mentioned after eqn (|4.7| ) the increase of 
the static energy due to the inclusion of the uj meson is no surprise although the number 
of degrees of freedom has been enlarged. The fact that at the minimum we always find 
A < 1 corroborates the repulsive character of the uj meson. Now we are finally able to 
mention that the test profiles, which we have been using all the time, correspond to the 
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parameterization ( p. 17] ) for m = 400MeV and A = 0.86, i.e. a configuration which seems 
to be very close to the stationary solution. Table 5.3 also contains the contribution of the 
uj meson to the energy as defined in eqn ( [5.16| ). The fact that this quantity is more sensi- 
ble to the chosen parametrization than the total energy indicates that a parameter ansatz 
for the soliton profiles is suitable to compute the minimal energy but some refinements 
may be necessary when studying other quantities, e.g. radii, magnetic moments, etc. . 
In figure 5.2 the profile functions which minimize E^^^^ in the parametrization ( |5.17| ) are 
plotted. In this figure the quark baryon density, b ocRHS of eqn (|5.6|) , is scaled such that 
the spatial integrals over u/m and b coincide. Obviously the uj profile is enhanced over 
b at large distances leading to a larger isoscalar radius in agreement with assertion (|1.3|). 



The fact that the presence of the u meson causes the soliton configuration to swell has 
also the consequences that the axial charge becomes larger and that the valence quark 
gets bound more strongly. These results were also obtained in the approach based on the 

We have furthermore considered fluctuations around the 



linear approximation ( 4.12 
profile 9 k 
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:jmm(r) lu ordcr to obtain an improved solution. The best we were able to achieve 
was a further decrease of the total energy by about 5MeV, i. e. the parametrizations ( |5.17| ) 



and ( 5.18 ) allow for a good approximation to exact minimum. 

We should also note that we have performed the analogous calculations using the 
regulator function /. In that case no local minimum was obtained. This provides a 



further indication why the approach of ref . |^T| did not lead to self-consistent solutions. 

Here we will not further discuss the results for various observables in this approach but 
rather make some remarks on the case when the isovector mesons p and ai are included. 
The previous formulas apply as well, however, the Hermitian part of the Hamiltonian has 
to be supplemented pi?[ . 



h 







(X ■ p + m(3 (cosO(r) + i'-^^T ■ f'sin0(r)) 

+-{a. X r) ■ TG{r) + -{ct ■ r){T ■ r)F{r) + -(cr ■ T)H{r). 

Zi Zi Zi 



(5.19) 



The radial functions G{r),F{r) and H{r) denote the vector (V^ = V^ ■ t/2) and axial 
vector {Af^ = A^ • t/2) meson profiles 



V-"" 

I 



eakihG{r) AI = 5iaH{r) + hraF{r) (z, a = 1, 2, 3). 



(5.20) 



As the number of fields involved increases the parametrization of the profile functions 
entering the Hermitian part of the Hamiltonian get more ambiguous. Motivated by the 
breathing mode description one might, in addition to ( ^.17 ), assume the ansdtze 



$;,(r) = A"<l>s.c.(Ar) for $ = F, G, H. 



(5.21) 



The choice n = would be analogous to ( |5.17| ) while ra = 1 is suggested by the asymptotic 
form of the (axial) vector meson fields in terms of the chiral angle. 

Two important changes in comparison to the previous studies have to be made in 
order to compute the total energy when the p and ai fields are included. First the mixing 
between the pions and the axial vector meson field leads to a different relation between 
the cut-off A and the pion decay constant /^r IjT^ 



f 

J IT 



Ncm^ 

47^2 




(5.22) 
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Table 5.4: The position (A™™) and the value (^E^^^^^) of the local minimum of the energy 
when the vector mesons are included. Also shown is the contribution of the u field (|5.16|) . 
The constituent quark mass m = 350MeV is assumed. 



n 


^min 


^Mmk (^^y) 


E^ (MeV) 


0.00 


0.850 


1662 


553 


0.25 


0.800 


1622 


490 


0.50 


0.745 


1565 


418 


0.75 


0.680 


1501 


340 


1.00 


0.630 


1433 


285 


1.25 


0.600 


1387 


253 


1.50 


0.603 


1376 


249 


1.75 


0.625 


1385 


259 


2.00 


0.655 


1403 


273 



This increases the value for A once a constituent quark mass m is given. Secondly the 



vector meson mass term becomes 24 



2^M°'(a 




G'ir) + ^F^ir) + F{r)H{r) + '^H^r) - 2ujl{r) 



(5.23) 



We are now enabled to search for a local minimum using the parametrization ( ^.21 ). In 
table 5.4 we present values A™™ of the breathing mode variable which lead to a local 
minima as well as the corresponding energies in the interval < n < 2. Here we use the 
constituent quark mass m = 350MeV and assume the imaginary part in the unregularized 
fromP. As a reminder we should mention that without the u meson this soliton energy is 



obtained to be lOlOMeV |]T9|. We observe that a local minimal in the energy surface is 
obtained for n ^ 1.5 and A ~ 0.58. Upon further increase of n the energy increases rather 
slowly. In the same way the value of the scaling parameter corresponding to the minimal 
energy, A™°, gets larger again. The minimum of the total energy is also characterized by 
the smallest contribution from the u meson. It should be noted that the minimal energy 
(1376MeV) is still 326MeV above the energy of the trivial B = 1 configuration where 
three (almost) free valence quark orbits are occupied. Again we observe the repulsive 
character of the u meson reflected by the fact that A™™ < 1. In contrast to the case 
with only the pseudoscalar fields this leads to a somewhat smaller binding of the valence 
quarks. We find e^^i/m = —0.29 to be compared with —0.38 in the model without the uj 
||19|| . This reduction seems to be due the diminishing profile functions of the (axial) vector 
mesons. A similar behavior has been obtained when extremizing the energy functional 
involving the definition (|4.13|) in the case that the imaginary has not been regularized 



24| . As this prescription is also adopted here we recognize at least qualitative similarities. 



The value for the axial charge of the nucleon as obtained from the integral p^ 



9A = - 



2ti 

^92 



drr 



H{r) + -F{r) 



(5.24) 



increases from 0.27 to 0.48 by the inclusion of the u meson. This is, of course, a correction 
into the right direction but nevertheless only about 1/3 of the empirical value extracted 



^When p and ai fields are included the regularized baryon number deviates strongly from unity |24 
in contrast to the model with pseudoscalars only. 
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from neutron /?-decay. It has been argued that 1/Nc corrections may further improve 
on the resuh for ga |3^; these arguments have, however, to be taken with some care in 
respect to the symmetries of the model 



As the main subject of the present paper is the investigation of the u meson in the 
context of the energy functional we will not go into detail on the nucleon properties. 
Especially because we only have a parametrical description of the stationary field con- 
figuration rather than a self-consistent solution. The nucleon properties are most likely 
more sensitive to the shape of the meson profiles and the quark wave-functions than the 
soliton energy. Nevertheless the construction of the self-consistent solution would be an 
interesting path to pursue in particular in order to gain information about the nucleon 
properties. In this context we also wish to remark that the uo profile associated with the 
minimal energy configuration (n = 1.5, A™™ = 0.598) does not turn out to be a smooth 
radial function but rather possesses small fluctuations off a monotonously decreasing 
function. We expect an exact solution to improve in this respect. 

6. Summary and Conclusions 

In this paper we have performed a thorough study of the analytical properties of 
a fermion determinant for the case that the a;-field is present in the static hedgehog 
configuration. This has recently been a matter of dispute in the context of the soliton in 
the NJL model with (axial) vector meson degrees of freedom. The fact that this model 
requires regularization makes mandatory the continuation forth and back from Minkowski 
to Euclidean space. This continuation represents an element of the transformations we 
have been considering in the complex z-plane defined by the Dirac Hamiltonian (p^.4D . As 



originally the fermion determinant is expressed in terms of the eigenvalues of this non- 
Hermitian operator we first have investigated the behavior of the associated eigenvalues in 



the 2;-plane. Thereby we have confirmed an earlier result |]2B| that for physically motivated 
field configurations these eigenvalues are smooth functions along the path connecting 
Euclidean and Minkowski spaces. Furthermore we have demonstrated that the Laurent 
series of such an eigenvalue or -function indeed reduces to a Taylor series. In this way 
the analyticity of eigenvalues and -functions is verified. Such an analytic behavior is 
not a priori clear because the eigenvalues are roots of a polynomial of large degree. We 
have, however, encountered a different kind of non-analyticity which is caused by a level 
crossing. Such crossing may appear if the cj-field is sufficiently pronounced. Thus, in 
order to attach physical significance to a given field configuration one has to check that 
no such level crossing takes place along the path from Euclidean to Minkowski spaces. 

In the next step we have then examined the behavior of the full fermion determinant 
in the complex 2;-plane. Formally analytical one-particle eigenvalues yield an analytical 
determinant in the unregularized form. However, regularization is unavoidable and in- 
troduces (at least in the proper-time scheme) the Hermitian conjugate of the Euclidean 

Dirac operator p^^. As a consequence the real and imaginary parts of the one-particle 
eigenvalues are treated differently when summing up for the fermion determinant. This 
represents the main reason for the fermion determinant to be non-analytical in z. As a 
matter of fact the corresponding Laurent series (|3.25|J3.2^) has vanishing radius of con- 



vergence. We therefore conclude that the analytical continuation of the Euclidean action 
functional to Minkowski space does not exist. This makes a derivation of the Minkowski 
energy functional by means of analytic continuation impossible. 

In order to nevertheless obtain a sensible definition of a Minkowski energy functional 
we have then formulated two consistency conditions to be satisfied by this functional. 
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The first one says that the energy changes by NcBuo^ under a constant shift (cUc) of the 
cu-field as a consequence of global gauge invariance {B is the baryon number.). Secondly, 
the current field identity for the baryon current imposes a normalization on the cu-field 
( |4.7| ). We have then examined the three treatments known in the literature with respect 
to these conditions. In this context only two approaches (by Schiiren et al. [^ ^ and 



the one by the present authors p5| , |2^ ) were shown to be consistent. The reason is that 
these two approaches are motivated by the Euclidean energy functional in contrast to the 
treatment by Watabe and Toki |20[. 



In the last section we have then investigated the expansion of the fermion determinant 

in powers of the w-field. In particular we have examined the role of p^. For this purpose 
we have considered the w-field as a perturbation to the Euclidean Dirac Hamiltonian and 
applied standard perturbation techniques. We have then demonstrated that at second 
order the bilocal kernel K, in eqn ( ^.7|) is not positive definite putting some doubts on the 
validity on an approach which is based on counting u^^ powers in the eigenvalues e^. The 
reason is that the corresponding regularization function / [e^, e^; A) does not smoothly 
assume the limit / — >■ as e^ — > ej,. A vanishing regularization function for identical 
arguments is required in order to satisfy the consistency conditions derived in section 4. 

As a matter of fact the functional trace of the operator PePe rather than the eigenvalues 
of the Dirac Hamiltonian is required in order to evaluate the fermion determinant in 
the proper-time regularization scheme. Stated otherwise, the proper-time prescription 
has to be imposed at the operator level. This point of view has already been adopted 



previously when performing the semi-classical quantization for two ||2^ and three flavors 



p7| . We have then observed that imposing the proper-time prescription at the operator 



level leads to a different expansion in terms of the tu-field for the action functional. From 
a mathematical point of view there are two ways for understanding the different results. 
First, in order to obtain the regularized form of the energy functional ( p.l5| ) from the 
regularized action (|1.9|) one has to make the inexact assumption that the Euclidean Dirac 
Hamiltonian and its complex conjugate can be diagonalized simultaneously. It is then 
also clear that both expansion schemes yield identical results when the perturbation of 
the Dirac Hamiltonian is static and Hermitian. We have verified this for an expansion 
up to second order. The second explanation relies on the fact that due to regularization 
the Hilbert space gets restricted and thus the rules for handling logarithms are no longer 
valid. Adopting this explanation it becomes also clear that in the limit A ^ oo (no 
regularization) the two expansion schemes are again identical. In regard of the fact that 



the derivation of the functional ( |3.11| ) involves a non-converging quantity, the results 



for the semi-classical quantization p7| , p7| , the smoothness of the regulator function / 



( |5.12| ) and the existence of a local minimum in the energy surface ( [5. 141 ) we consider 



the imposition of the proper-time regularization prescription at the operator level as the 
correct one. This conclusion finally prohibits an expansion of the eigenvalues e^ in order 
to obtain the Minkowski energy functional in the presence of the uj meson. We have also 

seen that a series based on the expansion of PePe converges rather quickly (at least at 
second order). 

This conclusion has enforced us to investigate the corresponding energy functional 
( 5.14|) in more detail. Employing a parametrical description of the chiral angle we have 



constructed the local minimum to this functional. Local minima exist for constituent 
quark masses m^400MeV. We have great confidence that the one-parameter ansdtze 
(I5.17J5.18 ) describe the exact solution quite well. For these solutions we have confirmed 
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the repulsive character of the cu-field. When other vector meson fields entering the Her- 
mitian part of the Dirac Hamiltonian are included as well, a parametrical description of 
the Dirac Hamiltonian is more ambiguous. Nevertheless we have demonstrated that for 
the constituent quark mass m = 350MeV a solution exists. In particular, the result that 
the valence quark level joins the Dirac vacuum, which was previously obtained in calcu- 
lations ignoring the u meson |]T^, is affirmed. This is again mainly due to the repulsive 



character of the isoscalar u. 

Let us finally add a note on similar soliton models which include the u meson. In 



extended Skyrme models [|11|, ^, |3J, ^ or the chiral quark model [^ the u meson 
appears only in terms of explicit polynomials which renders the energy functional an- 
alytic in the complex ^-plane. More importantly these models do actually not require 
regularization abandoning the need for continuation forth and back from Minkowski to 
Euclidean spaces. 

In order to point out an interesting task to pursue in future we should mention that 
we have not yet solved the complete set of stationary conditions in the presence of all 
meson fields. In the present work we have only made plausible that such a solution exists. 



However, a calculation aimed at the extremization of the functional ( 5.14 ) with all meson 
profiles included, which are compatible with the grand spin and parity symmetries of 
the hedgehog ansatz, is highly desirable because we have (for the first time) available 



a reasonable approximation (|5.14 ) to the energy functional of the NJL soliton model 



including all relevant meson fields. This should represent a well suited model for the 
description of static nucleon properties. 
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Appendix: Regularization for Quadratic Terms 

In this expansion we wish to describe the expansion of the vacuum part of the Eu- 
clidean energy functional in terms of the isoscalar vector cj-meson. The starting point is 
represented by the Euclidean Dirac operator 

i(3pE = —dr — h = —dr — ho, — iuj^ (A.l) 

with the Dirac Hamiltonian being decomposed into Hermitian {ho) and anti-Hermitian 
{iuja) pieces, see also eqn (^.1]). In order to extract the vacuum part we will only consider 
the limit of large Euclidean times T — > oo. We actually need to expand 

Ae = Tr log {-dr -he- iuj^) (A.2) 

in powers of the a;4 field: Ae = Ae + Ae + A^e + • • •• Here the superscripts label the 
order in which the u) field appears. The term linear in 004, is straightforwardly obtained 
to become purely imaginary 

Al = -zTt [{dr + he)-^ 004] . (A.3) 
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In the limit of large Euclidean times (T -^ oo) the temporal part of the trace is replaced 
by a spectral integral while the trace over the remaining degrees of freedom is carried out 
using the eigenstates of Hq {cf. section 3). We would like to remark that the corresponding 
eigenvalues e° are real. Then one easily finds 



,. J Z7T ,, J- 



°° dv 



-oo27rt;2 + (e 



0^2 



(z/|u;4|z/).(A.4) 



Again one may introduce a finite proper-time cut-off A by expressing the denominator 
in eqn ( A.4 ) in terms of a parameter integral 



/oo fj-i; roc |- , , -, 

-/ dse^v[~sU + {elf)]{v\u,\v). 
^ -oo zvr J1/A2 L \ yj 

The spectral integral may now be performed which results in 



(A.5) 



^^ = YiV^r^sgn(6°)erfc(^ 



(z/|a;4|z/). 



(A.6) 



At this point is important to remark that this result can also be gained by the application 
of standard perturbation techniques to the expression for the Euclidean energy functional 
obtained in section 3. In a first step one expands the functional ( |3.24]) ioi z = i and 
z* = —i in terms of the one-particle eigenenergies. Secondly, standard perturbation 
theory is employed to make the u^ dependence explicit as described in eqns { f).l\ j and 
( |5.2| ). Carrying out this procedure up to linear order in 004^ exactly reproduces eqn ( [A .61 ). 
The term of quadratic order el will be considered later. 

Next we turn to the evaluation of A%. As the terms of even powers in u^ contribute to 
the real part of the fermion determinant we may start the definition of the proper-time 
regularization (|1.9| ) 



A 



R 



-It. 

2 



dS / Tft Tftt 

— exp -sPePe 

1/A2 S \ 

4 



(A.7) 



In the presence of the uj meson the operator PePe is presented in eqn (|5.9| ). We then apply 
techniques which have previously been developed for the description of meson fluctuation 



off the solitonpQf to this expression. As the oj field is assumed to be static the expansion 
is also similar to the 1/A^c expansion performed to extract the moment of inertia in the 
cranking approach [^. The non- vanishing contributions at quadratic order in lo^ are 
found to be 



A^ 



R 



1 r°° 

-Tr / ds exp 

2 A/A2 ^ 



— Tr 



{H 





1-x 



dl 



UJa 



(A. 



xexp 



' dss dx dy exp —sx (/i|, — df) {2iuj4dr + i [he,uj4]} 

1/A2 Jo Jo L V n 

-s(l -x-y){hl- dfj {2iuJidr + i [he, wj} exp -sy {h% - dl 



One of the Feynman parameter integrals can be carried out trivially due to the cyclic 
properties of the trace. Furthermore the temporal part of the trace turns into a spectral 
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integral for T ^ oo. 






-Ttr / 
2 Jl/A- 



as / — exp — s I aIq + 



t; 



t^. 



(A.9) 



— -Ttr / dss I — I da exp 



1/A2 



27r 







—salkQ + v"^) {2UJ4V + i[hQ,uj4\} 



xexp 



s(l — a) (h^ + v^) {2iuj/^v + i [ho,, 074]} . 



These integrals are Gaussian and may readily be carried out. The remaining trace (tr) is 
computed using the eigenstates of /le- These manipulations then result in 



Al 



Nc^^ [^ ^ ds _,(,0)2, I I , 
00 ds 2 + s(e?-e0^2 



-^ty / y ^ 



(A.IO) 



Nc 



UfJ, 



1/A2 y/Ans^ 



(^y-(J>)i (e-"'°''-e-'"?'')(H^.I^>(A.|^.k>. 



Making use of closure as well as of the invariance of the second term under e° ^->- e° this 
expression may be rewritten as 






X 



ds 



1/A2 a/4 



TTS 



■.{y\u)^\li){\i\u)^v) 



,0e-s(e",)2^^0g-.(el^ ^g-s(.0) 



°^' P-«(^S)^ 






^ (e°P-(e°P 



(A.ll) 



from which / (e^, e^; A) in eqn ( [5.12|) can be extracted after performing an integration 
by parts in the last term. Next the question arises whether the application of standard 
perturbation techniques yields the same result. As already discussed in the main part of 
this paper it does not! In order to see how this comes about we start from the expression 
of the real part of the action ( 3.14| ) and perform the spectral integration 



A 



N, 



c 



R 



ty: 



°° ds 



2 V-^i/A2 V47rs3 
Using the expansion ( ^.11) for the energy eigenvalues (e 



exp(-s(e^)2 



^ = 6° + e2 



Al = Ar^T5:/^^^6°e^exp(-.(6^y 



and substituting the expression ([5.2|) gives 



Al 



{u\uj4\fi){fi\uj^\u) 



Nc^ V- ^Sn(e°)erfc (|f |) - sgn(eO)erfc ( 









(A.12) 



(A.13) 



(A.14) 



For the energy functional this leads to the expression / (e^, ty\ A) given in eqn (^.5|). 



31 



References 

[1] For an extensive compilation of articles on soliton models for baryons see, e.g. 
"Baryons as Skyrme Solitons", G. Holzwarth (Ed.), World Scientific Publ. Comp., 
Singapore 1993. 

[2] S. Brodsky, J. Ellis and M. Karliner, Phys. Lett. B206, 309 (1988); 
J. Ellis and M. Karliner, Phys. Lett. B313, 131 (1993); 
J. Schechter, A. Subbaraman and H. Weigel, Phys. Rev. D48, 339 (1993). 

[3] G. 't Hooft, Nucl. Phys. B72 (1974) 461. 

[4] E. Witten, Nucl. Phys. B160 (1979) 57. 

[5] T. H. R. Skyrme, Proc. R. Soc. A260 (1661) 127. 

[6] G. S. Adkins, C. R. Nappi and E. Witten, Nucl. Phys. B228 (1983) 552. 

[7] E. Witten, Nucl. Phys. B223 (1983) 433, 442. 

[8] O. Kaymakcalan, S. Rajeev and J. Schechter, Phys. Rev. D31 (1985) 1109. 

[9] J. Schechter, Phys. Rev. D34 (1986) 868. 

[10] M. Bando, T. Kugo and K. Yamawaki, Phys. Rep. 164 (1987) 217. 

[11] Ulf-G. MeiBner, N. Kaiser and W. Weise, Nucl. Phys. A466 (1987) 685 (1987); 
Ulf-G. MeiBner, Phys. Rep. 161 (1988) 213; 

P. Jain, R. Johnson, Ulf-G. MeiBner, N. W. Park and J. Schechter, Phys. Rev. D37 
(1988) 3252; 
Ulf-G. MeiBner, N. Kaiser, H. Weigel and J. Schechter, Phys. Rev. D39 (1989) 1956. 

[12] J. J. Sakurai, "Currents and Mesons", University of Chicago Press, Chicago 1969. 

[13] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 246. 

[14] For a recent review see e.g. T. Hatsuda and T. Kunihiro, "QCD Phenomenology 
based on a Chiral Effective Lagrangian". Phys. Rep., to be published. 

[15] D. Ebert and H. Reinhardt, Nucl. Phys. B271 (1986) 188. 

[16] J. Schwinger, Phys. Rev. 82 (1951) 664. 

[17] H. Reinhardt and R. Wiinsch, Phys. Lett. B215 (1988) 577; Phys. Lett. B230 (1989) 
93; 

T. Meissner, F. Griimmer and K. Goeke, Phys. Lett. B227 (1989) 296; 
R. Alkofer, Phys. Lett. B236 (1990) 310. 

[18] R. Alkofer and H. Reinhardt, Phys. Lett. B244 (1991) 461. 

[19] R. Alkofer, H. Reinhardt, H. Weigel and U. Ziickert, Phys. Rev. Lett. 69 (1992) 
1874; 
F. Doring, E. Ruiz-Arriola and K. Goeke, Z. Phys. A344 (1992) 159. 



32 



[20] T. Watabe and H. Toki, Prog. Theor. Phys. 87 (1992) 651. 

[21] C. Schiiren, E. Ruiz-Arriola and K. Goeke, Phys. Lett. B287 (1992) 283. 

[22] C. Schiiren, F. Doring, E. Ruiz-Arriola and K. Goeke, Nucl. Phys. A565 (1993) 687. 

[23] R. Alkofer, H. Reinhardt, H. Weigel and U. Zuckert, Phys. Lett. B298 (1993) 132. 

[24] U. Zuckert, R. Alkofer, H. Reinhardt and H. Weigel, Nucl. Phys. A570 (1994) 445. 

[25] S. Kahana and G. Ripka, Nucl. Phys. A429 (1984) 462. 

[26] C. Schiiren, Contribution at the ETC* workshop on The Quark Structure of Baryons 
in Trento, October 1993, unpublished. 

[27] H. Reinhardt, Nucl. Phys. A503 (1989) 825. 

[28] R. Alkofer, H. Reinhardt and H. Weigel, "Baryons as Chiral Solitons in the Nambu- 
Jona-Lasinio Model" , in preparation. 

[29] H. Weigel, H. Reinhardt and R. Alkofer, Phys. Lett. B313 (1993) 377. 

[30] G. S. Adkins and C. R. Nappi, Phys. Lett. B137 (1984) 251. 

[31] N. S. Manton and P. J. Ruback, Phys. Lett. B181 (1986) 137. 

[32] M. Wakamatsu and T. Watabe, Phys. Lett. B312 (1993) 184; 

A. Blotz, M. Praszalowicz and K. Goeke, Phys. Lett. B317 (1993) 195. 

[33] R. Alkofer and H. Weigel, Phys. Lett. B319 (1993) 1. 

[34] H. Reinhardt and D. Ebert, Phys. Lett. B173 (1986) 459. 

[35] B. Schwesinger, H. Weigel, G. Holzwarth and A. Hayashi, Phys. Rep. 173 (1989) 
173. 

[36] T. Cohen and M. Banerjee, Phys. Lett. B230 (1989) 129. 

[37] H. Weigel, R. Alkofer and H. Reinhardt, Nucl. Phys. B387 (1992) 638; 
A. Blotz et al., Nucl. Phys. A555 (1993) 765. 



33 



Figure captions 



Figure 2.1 

Left: The behavior of the valence quarks' eigenenergy in the complex plane z = Se^'^. 
The starting point </? = is indicated. From there it continues anti-clockwise as (f 
increases. Right: The same for the state with the smallest (positive) real part of the 
energy eigenvalue in the G^ = 2"'" channel. 



Figure 2.2 

The behavior of the valence quarks' eigenfunction in the complex plane in arbitrary units. 
ip labels the phase of the point in the complex plane z = 5e^'^ . Left: At r = 0, only the 
upper component is non-zero. Right: At an arbitrary intermediate point for the case 
6=1. 



Figure 3.1 

The behavior of the energy ii^jet (3.24) in the complex plane for the parametrization 
z = 5e^'^ . The starting point v? = is indicated. From there the energy continues 
anti-clockwise. 



Figure 5.1 

Left: The total energy as a function of the scaling parameter A defined in eqn (5.17) 
Right: The parametrical dependence of the energy on the parameter Rm as defined in 
eqn (5.18). 



Figure 5.2 

The profile functions which minimize the energy functional (5.14) under variation of the 

scaling parameter A. The quark baryonic density h ~ gl g {RHS of eqn (5.6)) is artificially 
scaled such that the spatial integrals over u/m and b coincide. 
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